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Abstract 

We calculate the specific heat for an interacting Fermi system near the fer- 
romagnetic phase transition using the Renormalization Group method. The 
temperature dependence of the specific heat present for dimension D = 3 a 
logarithmic dependence which shows that the fermionic excitations reaches a 
non - Fermi behavior. The result is in good agreement with the experimental 
data obtained recently for Ni x Pd\^ x alloys. 
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I. INTRODUCTION 



The occurrence of the non - Fermi behavior in the heavy - fermion system has been 
associated with the antiferromagnetic quantum phase transition (AQPT) [|T]-[|] and is one of 
the most important subject concerning the non - Fermi behavior, introduced for the expla- 
nation of the normal state of cuprate superconductors. Recently the non - Fermi behavior 
was experimentally discovered in the systems as Th\_ x XJ x CuiSi<i [[§] and the Ni x Pd\_ x alloy 
||, if x < x c , x c being a critical concentration. This behavior has been also considered 
as a quantum phase transition (QPT), but the systems present a ferromagnetic quantum 
phase transition (FQPT). In this paper we will calculate the heat capacity for a model 
which describes the interaction between itinerant - electron and magnetic moments. The 
Renormalization - Group method (RG) used in previous papers |]-|| will be adopted to 
describe the QPT and the free energy calculated by RG will be used to show the existence 
of the logarithmic term (T In T) which represent the non - Fermi behavior of the fermionic 
excitations. 

The paper is organized as follows. In Sec. II we present the model and give the dynamic 
susceptibility of the system. Sec. Ill contain the RG scaling equations and in Sec. IV we 
calculate the heat capacity. The relevance of the results compared with the other models is 
discussed in Sec. V. 



II. MODEL 

The transition to the ferromagnetism due to the magnetic moments has been discussed 
by Doniach and Wohlfart [[J for the palladium ion. The basic idea of the model is the 
polarization of the itinerant - electrons by the magnetic impurity and the susceptibility was 
calculated as: 

/ \ xo(q,w) 
x(q, «) = : r. 2J2R ,i ; : (1) 
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where x (q, cu) is the susceptibility of electrons given by 



xo(q,w) 
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xo(q,w) = —J2 



Tiki — n-k+ql 



(2) 
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A is the gap introduced by the electron - electron interaction J, and the electron - impurity 
interaction J, and was calculated as: 



A = IR + 2JR' 



(3) 



where 



R = jy S( n *i - n *T) 



and 



(4) 



R' = -C <S Z > 
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and C is the concentration of the magnetic impurities with the spin S. 
Using the expansion: 
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A, B, C, Di (i = 1,2) being constants we get for susceptibility expressed by Eq. ([I]) the 
form: 



x(q,^) 
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where 5 is the concentration dependent parameter, given by: 



s = x(o,oy 1 -F 



(8) 



F = I + 2J 



R 



(9) 



For A = Eq. (J7p approximate the dynamic susceptibility of the excitations in the itinerant 
ferromagnet (see Ref. ||[7|)- 
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III. SCALING EQUATIONS 



The fluctuations of the magnetization in the critical region of the ferromagnetic state 
are described by the action: 



where 



C _ c( 2 ) , c( 4 ) 
Beff ~ Jeff + b eff 



(10) 



(11) 



^ i=l qi \i=l J 
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where 



E = ^E 



(13) 



nJ (2vr) d 

and g» = (qi,u ni ), u n i = 2nnT being the bosonic frequency and u the interaction, which 
satisfies u > 0. Next we consider the scaling properties of this action. When we take the 
scale transformation q' = ql and u' = uil z , the Gaussian term will be transformed as: 

1 



c(2) _ * V- 
D e// - 2 
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The fourth - order term, given by Eq. ( |T2"D is scaled as: 

1=1 9 ' 4 \i=l / 
Therefore, the scaling equations should be: 

dT{l) 



(15) 



dA(Z) 
dl 



zT{l) 



(4-z)A(i) 



(16) 



(17) 



dl 

dS(l) 



,n:i ' ] -(z-3)T(l) (18) 



dl 



25(l) + 2u(l)(n + 2)fx (19) 



Ml = [4 - (d + z)\u{l) - (n + 8)u 2 (l)f 2 (20) 

The Eq.(^) shows that the quadratic term is scaled to zero for d + z > 4. On the other 
hand the existence of the ferromagnetic state is conditioned by the existence of the gap A 
and from Eq. (|17|) we may conclude that z < 4. However, from fll8|) we see that for z > 3 
the quantum fluctuations are important and for z = 3, T is dangerous irrelevant. These 
considerations leads to the conclusion that for d = 3 the relation d + z > 4 and u indeed 
scales to zero. 

If only the Gaussian point of action 5^ is considered, the partition function is 



Z = Trexp(-/5F) = J D<j>*D<j>exp{S) 



where the free energy F is given by: 

F = V f A /lf rk ^n B (u;) tan- (21) 

nsioj) being the Bose function. Performing the same scaling for q and u we obtain the 
equation for F: 

dF(l) 



dl 



(d + z)F(l) + f 3 (T(l),A(l),S(l)) (22) 



Following the method developed in Ref. |7||| we will calculate the free energy and the 
specific heat for the two regimes which exist when we stop the scaling. These regimes are the 
quantum regime (T(Z) 1) and the classical regime (T(Z) ^> 1) and we expect an important 
enhancement of the behavior in the specific heat due to the quantum effects. 
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IV. SPECIFIC HEAT 



The specific heat C v = —Td 2 F/dT 2 will be calculated from the general solution of the 

Eq. 

FIT) = f dxe-^ x f 3 (Te zx ) (23) 
Jo 

following the method used in [§],[7]]. 

The effect of the fluctuations is important in the quantum domain (T(/) <C 1) and 
classical domain (T(Z) ^> 1). These two domains have been studied using the flow equations 
(see Ref. 0) and the basic point is to stop the renormalization procedure when the system 
is driven from the critical region in the region where we can perform a standard perturbative 
calculation. We introduce I = I* which near the critical point satisfies T(l*) <C 1 for quantum 
domain and T(l*) ^> 1 for the classical domain. This domains are separated by a matching 
value l m defined by: 

T(l m ) = 1 (24) 

calculated from Eq. (16) as: 

l m = \ In ^ (25) 
In this approximation we calculate the free energy from Eq. ( p3l) as: 

F(T) = / dxe- {d+z)x h (Te zx ) + / dxe~ {d+z)x f 3 (Te zx ) (26) 

J Jl/zlnl/T 

Following Ref. we approximate f 3 {T) — f 3 (0) ~ AT 2 in the low temperature domain 
and fs(T) ~ BT in the high temperature domain. The general form of the free energy given 
by (H) is: 

F(T) = F{F 2 In 1 - F 2 T (27) 
and the specific heat obtained from Eq. ( p7|) has the form: 
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C v = 7 oT + 7l TlnT (28) 

In this equation the first term represent the effect of critical fluctuations in the classical 
domain and the second term describes the effect given by the quantum fluctuations in the 
proximity of the quantum critical point, where the system becomes non - Fermi in agreement 
also with the resistivity temperature dependence 



V. DISCUSSIONS 

We will discuss the results obtained in this paper in connection with the other approaches 
as well as to the experimental results. 

First we mention that recently the problem of QPT in a ferromagnet was extensively 



studied |]TT|-|T4|| for a pure as well as for disordered ferromagnet. The main idea of these 
papers is the occurrence of an effective long - range interaction between the order parameter 
fluctuations. In fact we mention that this conclusion is valid only at T = and the non - 
analyticity is done to the fact that in JTT|-p^] the authors integrated out all the fermionic 
degrees of freedom [|T5[]. However, in this case the coupling to the electron - hole continuum 
is not any more possible, and the existent damping term is difficult to be explained in pure 
ferromagnet, at least due to this effect. 

An important point is that such a non - analytic susceptibility does not respect the 
conservation law for the magnetization in an itinerant - electron ferromagnet. The difficulties 
due to the application of the RG method as well as the conjecture that for T^O the model 
becomes massive by the simplest substitution — > (q + T) d ~ x determined us to use 
the Hertz - Millis (HM) approach in description of the ferromagnetic quantum phase 
transition. 

The formal inconvenience which appears by taking into consideration the occurrence of 
many energy scales seems to be normal, but also natural if we are crossing from quantum - 
to classical regime. 

The occurrence of the non - Fermi behavior in this model is also a good reason to use 



the HM version of RG for this problem. The only problem which remain open is the value 
of z in this model, but we expect that the value from the disordered ferromagnet to be more 
appropriate. 

Finally, we mention that we obtained a good agreement with the experimental data JT]|| 
for the specific heat and we also need an accurate calculation for the magnetic susceptibility 
as a function of temperature. A preliminary calculation showed a strong deviation from 
normal Pauli behavior, but also from the Curie - like behavior as is expected from self - 
consistent fluctuation theory ||10|| . The occurrence of a In T - term predicted by experimental 



results in AQPT is due to the non - Fermi behavior. The accurate analysis of this problem 
for FQPT is in progress. 
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